Mon. Not. R. Astron. Soc. 000, ITHT31 (20081 Printed 1 December 2009 (MN WT$i style file v2.2) 



Sizes and ages of SDSS ellipticals: Comparison with 
hierarchical galaxy formation models 
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ABSTRACT 

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that 
the shape, normalization, and dispersion around the mean size-stellar mass relation 
is the same for young and old systems, provided the stellar mass is greater than 
3 x 10 10 M Q . This is difficult to reproduce in pure passive evolution models, which 
generically predict older galaxies to be much more compact than younger ones of the 
same stellar mass. However, this aspect of our measurements is well reproduced by 
hierarchical models of galaxy formation. Whereas the models predict more compact 
galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the 
more massive objects, resulting in a flat size-age relation at the present time. At lower 
masses, the models predict that mergers are less frequent, so that the expected anti- 
correlation between age and size is not completely erased. This is in good agreement 
with our data: below 3 x I0 10 M©, the effective radius R e is a factor of ~ 2 lower for 
older galaxies. These successes of the models are offset by the fact that the predicted 
sizes have other serious problems, which we discuss. 
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1 INTRODUCTION 

According to the standard cosmological paradigm of 
structure formation and evolution, dark matter (DM) 
halos have grown hierarchically, through the merg- 
ing together of smaller units into ever l arger systems 
dPress fc Schechterl [l97i : iBond et all Il99ll : lLacev fe CoTel 

galaxies form 
of DM halos 



and active galactic nuclei (AGN) 


(e.g., 


Benson et al. 


2003; Granato et al. 2004, 2006: 


Di Matteo et al. 


2005; 


iBower et al.ll2006l: ICattaneo et al.ll2006l; ICroton et al.1 


2006; 


De Lucia et al. 


2006; 


Menci et al. 




20061; 


Monaco et al. 



1 19931 ; ISheth et al.l l200l"h . In this scenario 
inside this hierarchically growing system 



l|Whitel Il979l ; IWhite fc Frenklll99lh . Sem i-analytical mod 



els (SAMs) of ga l axy formation (e.g. , I Cole et al.l 



2004 



Benson et aj1l2003l; iGranato et~ail 12004 iMenci et al 
Bower et all 120061; ICattaneo et al.l |2006|; ICroton et al.1 



2000 



20061; iDe Lucia et al.1 l200rj; IKhochfar fc Silkl l2006al lbl: 



Hopkins et al. 20061 ; Monaco et al.l 12007 ~) have now been 



2007). While this feedback significantly reduces the amount 
of cooling baryons in the host halos, it also seems to 
be a promising tool to account for the AGN luminosity 
functions a nd mean trends in the stellar-halo mass rela- 
tions fe.g.. IWvithe fc Loebll2003l; IScannapieco fc Qhl I2004J: 



Cirasuolo et al.ll2005l; | Di Matteo et alj |200." 
20051; IVittorini et alT~ 



120051; 



Sazo nov et al.l 
200a 



Hopkins et al. 

Lapi et al.1 120061; IShankar et al.l 120061 ; iBenson et all 120071: 



20081 ; IShankar et al.l l2008b) 



Fontanot et al.1 120071; iMalbon et al.l 120071 ; iMarulli et all 



able to reproduce several properties of the local galaxy 
population. In particular, the local galaxy luminosity 
function has been successfully matched both at the high 
and low luminosity ends owing to the implementation 
of models for the feedback from both stellar evolution 
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In addition to the luminosities, masses, and metal- 
licities of galaxies, their sizes, measured at low and high 
redshift, pro vide strong constraints on galaxy formation 
models fe. g., ICirasuolo et al.1 120051: IKhochfar fc Silkl l2006a| : 



Almeida et al.l 120071: iFan et al.ll200Sl; iGonzalez et al.ll2009l 
Hopkins et al"f l2009bl~lBernardil 120091; Ivan der Wei et al.l 



2009). Recent observations show t hat galaxies of a given stel- 



lar mass are more compact fe.g.. lTruiillo et al.ll200rj. 120071 ; 
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Cimatti et al. 20081; Buitrago et al.l 20081; Chapman et al 



2008] ; iFranx et al.l |200S|; |Saracco et alj|2009l; iTacconi et al 



2008 



2008; van der Wcl ct al. 2008; van Dokkum ct al. 
Younger et al.l 120081 ; iDamianov et al.l 120091 ; I Williams et al.l 



2009), and, at least some, have higher velocity dispersi ons 
(e.g.. ICenarro fc Truiilloll2009l : Ivan Dokkum et al.ll2009l ) at 
redshifts greater than unity; this is also tru e of BCGs, even 
for small lookback times l|Bernardil 120091 ). [ Saracco et al.l 
(2009) also showed that at z ~ 1.4, while older galaxies 
at fixed stellar mass tend to lie a factor of ~ 2 — 3 below 
the size-s tellar mass relati on characterizing local early- type 
galaxies l|Shen et all 120031 ) . younger galaxies are consistent 
with it. After anal yzing a sample of 12 v ery massive galaxies 
at similar redshift, iMancini et alJ (2009) suggested a down- 
sizing scenario in sizes, with the most massive galaxies ap- 
proaching the local size-ma s s relat ion earlier then less mas- 
sive ones. ICappellari et al] (|2009l ) also discussed that two 
massive galaxies with individual spectra at 1.4 < z < 2 are 
similar to local counterparts, and another 7 galaxies with 
velocity dispersion from staked spectrum, are consistent to 
the most dense local galaxies of the same mass. 

Small sizes (and, at fixed mass, higher velocity disper- 
sions) at high redshifts are not unexpected: they result if 
galaxies at higher redshifts formed through more gas-rich 
and dissipative mergers (e.g.. iRobertson et al.ll2006h . More- 
over, the gas fractions and t he overall density of the uni- 
verse decrease with time (e.g.. lHopkins et alj2009bl . and ref- 
erences therein) implying, on average, less dense remnants 
at later times. Therefore if galaxies continuously form at 
different epochs, it is expected that in deep observation at 
2 > 1.5 — 2, galaxies at fixed stellar mass might follow a dis- 
tribution of sizes, with older galaxies born from gas-richer 
events, being more compact. 

However, the problem is to explain how all galaxies, 
irrespective of their age, evolve in time on the same loca l 
size-mass relation. A s show n by Ivan Dokkum et"al] (2008), 
IShankar fc Bernardil l|2009h . and further discussed here, a 
simple pure monolithic collapse followed by strictly passive 
evolution is not satisfactory, as this does not explain why the 
extremely small sizes and high densities of massive galax- 
ies at high-2 have no local counterparts in SDSS (also see 
iTruiillo et al]|2009l ). 

Shankar & Bernardil l|2009h used a sample of about 
48,000 early-type galaxie s from the Sloan Digital Sky Survey 
(SDSS; lYork et al] |2000) to show that older galaxies have 
smaller half-light radii R e and larger velocity dispersions a 
than younger ones of the same stellar mass Mstar- Specif- 
ically, they found that, using the age-corrected luminosity 
L£ orr as a proxy for Mstar, galaxies with age ~ 11 Gyrs 
below L£ orr ~ 10 11 Lq have R e smaller by 40% and a larger 
by 25%, compared to galaxies that are 4 Gyrs younger. The 
sizes and velocity dispersions of more luminous galaxies vary 
by less than 15%, whatever their age, implying a signifi- 
cant break in the R e -L^° TT and a-L^ 011 relations at high 
L£ orr . The galaxies in their sample have been carefully se- 
lected to be early-type galaxies, and they are all charac- 
terized by much larger sizes and lower densities than their 
higher redshift massive counterparts. Also, the differences 
between the sizes o f old and young galaxies reported by 
IShankar fc Bernardil l|2009h . are far less than expected from 
a simple monolithic model evolution with R e <x (I + 2) -1 at 



fixed stellar mass, which would result if the galaxy density 
is proportional to the density of the universe. 

Our goal in what follows is to compare these findings 
data with hierarchical models of galaxy formation and evo- 
lution to understand how well they can reproduce the data. 
In Section [2] we first revisit the main observational results 
regarding the size-mass relation in local ellipticals, showing 
that the very massive galaxies of any age follow similar size- 
stellar mass relations in shape, normalization, and disper- 
sion around the mean. Section[?]co mpares our measure ments 
with the publicly available SAM of iBower et"aH l|2006l . B06 
hereafter), which has successfully reproduced several statis- 
tical properties of galaxies. The B06 model, which is based 
on t he Millennium Simulation of the dark matter distribu- 
tion jSpringell 12005). provides the sizes of spheroids at any 
epoch, thus enabling a direct comparison between the pre- 
dicted and observed size-mass evolution. We use the SAMs 
to discuss how the number and type of mergers (dry or wet) 
scale with final stellar mass. In particular, we show that dry 
mergers, defined to be mergers between gas-poor progeni- 
tors, might be good candidates for erasing the effects of a 
monolithic collapse and producing a rather flat size-age rela- 
tion, similar to what observed in the data. Our conclusions 
are in § [5] where we also discuss some serious failures of the 
models. 

Throughout this paper we adopt the cosmologi- 
cal parameters Q m — 0.30, £7a = 0.70, and h = 
#o/100 km s^Mpc" 1 =0.7. 



2 DATA AND MEASUREMENTS 
2.1 The sample 

We use the SPSS - based sample of early-type galaxies from 
iHvde fc Bernardil ((2009). This sample was constructed in 
a way to maximize the contribution of ellipt ical galaxies as 
briefl y described here (see the recent work bv lBernardi et akl 
2009 for further details). The galaxies in this sample are 
very well described by a deVaucouler profile in both the g 
and r bands (fracDev—1), have "early-type" spectrum, with 
eClass< (to minimize later-type contamination), and have 
an additional cut in axis ratio of b/a > 0.6. These galaxies 
were also selected to have velocity dispersion a higher than 
60 kms -1 , close to the dispersion limit of the SDSS spec- 
trograph, and less than 400 kms -1 , to avoid contamination 
from double/multiple superpositions (Bernardi et al. 2006, 
2008). The resulting sample, which contains about 47,300 
early- type galaxies, is distributed within the redshift range 
0.013 < 2 < 0.3, which corresponds to a maximum lookback 
time of 3.5 Gyr. Based on the re cent, detailed analysis per- 
formed by iBernardi et alJ (120091 ) , this sample has minimal 
contamination by disky SO galaxies, which makes it an ideal 
catalog to compare with galaxy evolution models of massive 
spheroids. The galaxies in the sample have apparent mag- 
nitudes 14.5 < m r < 17.5 (based on deVaucouleur fits to 
the surface brightness profiles). (The SDSS photometric pa- 
rameters for these objects have been corrected for known 
sky subtraction problems which affect bright objects.) Es- 
timated stellar mass es and ages for these objects are from 
iGallazzi et all (|2005l ). These are based on running a like- 
lihood analysis of the spectra that returns a mass-to-light 
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Figure 1. Size-stellar mass relation for the galaxies in our sample with ages as labeled in the upper left corner of each panel. The solid 
lines bracket the regions containing 30%, 65%, and 95% of the full sample, respectively. The orange, green, and blue coloured regions 
show the corresponding distributions when we restrict the sample to narrow bins in (luminosity-weighted) age. Older galaxies shift to 
higher R e , and MgTAR, but are not offset from the relation defined by the full sample, whereas the youngest objects tend to be offset 
towards larger R e and smaller a. 



ratio M S tar/L z (defined for a IChabrierl [2003 initial stel- 
lar mass function), which is in turn converted to a stellar 
mas s using the SPSS pet rosian z-band restframe magnitude 
fsee lBernardi et al 1 l2009l for a comprehensive discussion and 
comparison of such stellar mass estimates with other meth- 
ods). Our m ain results do not cha nge if we use the ages 
published by iBernardi et al. ( 2006 ) , which were computed 
by fitting the iThomas et al.l (|2005l ) a-enhanced models to 



the Lick index abso rption features. The age estimates of 
Ijimenez et al.l (|2007l ). derived from single ste llar population 
spect ral fitting, using the MOPED algorithm H eavens et al.l 
|2000| ). also yield similar results. Also, a possible bias might 
be introduced by the youngest galaxies in the sample with 
ages < 5Gyr (the young age could be due to weak star for- 
mation which makes the galaxy to appear relatively young) , 



comparable to the lookback time of the sample and close 
to the average error in age estimates. However, we have 
checked that cutting out these galaxies and repeating the 
analysis does not minimally alter the overall conclusions of 
the paper. 

Finally, we note that the age and stellar mass estimates 
come from the same algorithm, so they have correlated er- 
rors. However, this does not bias the results which follo w 
(for details, see IBerna rdi 2009; Sha nkar fc Bernardj|2009h . 

In what follows we wish to study the bulk of the early- 
type population. However, e.g.. iBernardil (|2009l ) has argued 
that Brightest Cluster Gala xies (BCGs) had unu s ual fo rma- 
tion histories, so, following IShankar fc Bernard (|2009l ). we 
remove them from our sample. Our final sample is composed 
of ~ 45, 700 galaxies. 
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Figure 2. Corresponding velocity dispersion-stellar mass relations with same format as Figure[JJ Similarly, older galaxies shift to higher 
crand MgxARi but are not offset from the relation defined by the full sample, whereas the youngest objects tend to be offset towards 
smaller <r. 



2.2 Results 

It is well known that more massive early-type galaxies have 
larger half-light radii. The dotted line (same in all four 
panels) in Figure [JJ shows this relation in our dataset, ob- 
tained from a linear fitting log(R e /5kpc) = (—0.63 ±0.01) + 
(0.53 ± 0.01) log(MsTA R /10 ' M R ), in very go od agreement 
with previous results from IShen et al.l l|200St ). Notice also 
that the spread in size at fixed mass is approximately ~ 0.6 
dex for masses below Mstar~ 3 x 1O 1O M0, decreasing at 
large mass. This behavior in the scatter as a funct ion of 
stellar mass was already note d by Shen et al. ||2003|). and 



we here further suggest, after IShankar fc Bernardil (|2009l ) 
and the discussion below, that this is mainly induced by the 
different gradients in the size-age relation in different stel- 
lar mass bins. The solid contours, the same in each panel, 
bracket the regions containing 30%, 65%, and 95% of the 
full sample, respectively. 

The orange, green, and blue coloured regions show the 



corresponding distributions when we restrict the sample to 
narrow bins in (luminosity-weighted) age. Notice that the 
older galaxies tend to populate the high-mass end of this 
relation (the shaded regions peak at higher Mstar in the 
bottom right panel than in the top left), i n qualitative agree- 
ment with the notio n of downsizing (e.g.. ICowie et al.lll996l ; 
iHeavens et al.l [2004) . To proceed in a more accurate anal- 
ysis, we have refitted the size-mass relation for the sub- 
samples of galaxies considered in Figure [JJ We found that 
while younger galaxies follow a significantly shallower corre- 
lation, R e oc Mstar°' 48±0 ' 01 , the oldest ones follow a steeper 
relation, with R e oc Mstar°' 65±0 ' 01 • These behaviors are 
mainly caused by the fact that below Mstar J; lO n M0, 
older galaxies tend to gradually have sizes that are up to a 
factor of about two smaller than younger ones of the same 
Mstar- Above 10 11 M©, the sizes are instead similar, what- 
ever the age. 

Figure [2] shows a similar analysis of the velocity dis- 
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Figure 3. Mean effective radius as a function of age for a given subsamples denned by stellar mass; histograms show the age distribution 
of each subsample. The bottom right panel shows the mean effective radii found by averaging over 100 simulations in which each time 
a subsample of 1000 galaxies is randomly drawn from the catalogue: a mere random selection of galaxies from the parent sample does 
not reproduce the generally flat R e — age relation at fixed stellar mass. The dotted line in each panel shows the mean size for the stellar 
mass bin considered, as predicted by the global -Re-AfgxAR relation (dotted line in Figure [TJ. 



persions a rather than the sizes. As for the sizes, the over- 
all scatter in a at fixed Mstar is independent of age, but 
decreases with increasing Mstar, lending further support 
to our above suggestion that what is actually driving the 
mass-dependent scatter in the size-m ass relation is corre- 
lated to differences in ages (see also IShankar fc Bernardil 
2009). And, to lowest order, the mean velocity dispersion- 
stellar mass relation, which we fit as log(cr/200 km/s) = 
(-0.288±0.014) + (0.27±0.01) log(A/ S TAR./10 10 M Q ), is the 
same in all age bins - the primary trend being that older 
galaxies are shifted to larger Mstar values. A more detailed 
analysis yields that, at the low mass end, young galaxies 
(which had slightly larger sizes) have slightly lower velocity 
dispersions, thus causing a slight steepening of the relation. 

Figures [3] and [4] show another view of the correlations 
between size, mass and age: the size-age relation for bins in 
Mstar- In agreement with the previous Figures and what 
discussed in § U (see also IShankar fc Bernardil 120091 ) , the 
lower mass galaxies (upper left panels) have sizes and ve- 
locity dispersions which show some significant trends with 
age, with older galaxies having smaller R e and larger a. At 



the top of each panel we also indicate, for better reference 
to the models discussed below, the redshift corresponding 
to the lookback time equal to the age. The dotted lines in 
the Figures show the mean values of size or velocity disper- 
sion predicted by the global i? e -MsTAR and ct-Mstar rela- 
tions (dotted lines in Figures [1] and [4} for the stellar mass 
defined in each panel. It is clear that while the lowest mass 
bins present a significant gradient with age (opposite sign for 
sizes and velocity dispersion) , this trend progressively disap- 
pears when moving to more massive galaxies, with almost 
no correlation between size and age for the most massive 
ones. 

Nevertheless, more massive objects are offset to larger 
sizes. As a result, when averaged over all masses, older ob- 
jects have larger sizes, except possibly for the oldest galaxies. 
The bottom right panels of Figures and U show the results 
of this exercise. When randomly selecting galaxies from the 
whole sample (the points are averages over 100 realizations 
with 1000 points each), we find an increasing size and veloc- 
ity dispersion with increasing stellar mass. Thus, the size-age 
relation is almost entirely due to the size-mass and age-mass 
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Figure 4. Same pattern as Figure [3] Mean velocity dispersion as a function of age for a given subsamples defined by stellar mass, with 
histograms showing the age distribution of each subsample. The bottom right panel shows the mean velocity dispersions from random 
selection of galaxies from the parent sample. The dotted line in each panel shows the mean velocity dispersion for the stellar mass bin 
considered, as predicted by the global cr- .Mstar relation (dotted line in Figure [2}. 



correlations (massive objects are older and larger). This is 
analogous to the color-magnitude relation being entirely due 
to the correlati ons between color an d luminosity with veloc- 
ity dispersion (Bcr nardi et al.|[2005t ). 

Similarly, at fixed Mstar, there is little correlation be- 
tween a and age, but, because more massive galaxies are 
offset to larger a, the result of averaging the a -age relation 
over all Mstar yields a strong t rend: the oldest gal axies 
have the highest values of a (e.g., Bcrnardi et al. 2005|). 

2.3 Some implications 

In a strict monolithic scenario for galaxy formation in which 
the age of the stars reliably traces the time the galaxy was 
assembled, older galaxies should have smaller sizes than 
younger galaxies of the same stellar mass. This is because 
objects which assembled earlier did so in a denser universe. If 
the galaxy density is proportional to the density of the uni- 
verse at the time of assembly, one expects R e cx (l + 2f or m) -1 
at fixed stellar mass. (One might expect the actual scaling 
to be even stronger, since dissipation in the baryonic clumps 
from which the stars form is expected to be more efficient 
when the density is higher, so objects which form at higher 
redshift should be even more compact and have even higher 
velocity dispersions.) However, the sizes of old and young 



galaxies in our sample are much weaker functions of age. 
This rules out models in which star formation and assembly 
are concurrent, and galaxies passively evolve thereafter. 

On the other hand, as reviewed in Section[T] a large por- 
tion of high redshift galaxies are observed to be more com- 
pact than their counterparts at low z. Therefore, some pro- 
cess must have altered the sizes and stellar masses of galax- 
ies since they formed. Our results suggest that, whatever 
the mechanism, it must be "fine-tuned" so as to yield the 
essentially fiat size-age and er-age relations (at fixed Mstar) 
shown in Figures [3] and [3] 



3 COMPARISON WITH GALAXY 
FORMATION MODELS 

In this Section we compare our observational data with the 
predictions of sem i-analytic galaxy formation models. See 
iParrv et al] (120081 ) and Shankar et al. (2009c) for a detailed 
description and comparison of such models, which follow 
the cosmological co-evolution of dark matter halos, subha- 
los, galaxies and supermassive black holes within the con- 
cordance ACDM cosmology. Briefly, these models track disk 
and bulge components of each object as it evolves. B06 also 
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Figure 5. Predicted sizes (upper panels) and densities (lower panels), as well as the dispersions around the mean for galaxies of the 
same stellar mass, as labeled in each column, but identified at different epochs, in the Bower et al. (2006) model. On average, galaxies 
of the same stellar mass are predicted to be significantly smaller and denser at early times. 



output the half mass radius of the spheroid component, so 
this is the model we use as a reference in what follows. 

We will be interested in the masses, (half-mass) sizes, 
and formation histories of the early-types galaxies in these 
models. We classify a galaxy as early-type if it has Mbui ge — 
Mtotai < 0.4, where Mb u i gc and M to tai are the predicted 
B-band magnitude of the bulge and of the whole galaxy, 
respectively. This cut selects galaxies with a bulge-to-total 
ratio B/T > 0.7, which is the minimum B/T c haracter- 
izing the galaxies in the iHvde fc Bernardil {2009) sample. 
Note that this cut in B/T preferentially selects spheroid- 
dominated galaxies, while a more common cut of B/T > 0.5, 
comparable to a cut in concentration of C r > 2.86 (e.g., 
IGonzalez etaHl2009h . allows for a large c ontamination o f 
disky galaxies (see the recent analysis by iBernardil [2009). 
The latter type of galaxies might have had quite different 
formation histories than the ones of interest to this paper. 
In fact, a large fraction of SO galaxies in the B06 model 
grow their bulge via disk instability. Also, as further dis- 
cussed in Shankar et al. (2009c), the mock galaxy sample 
considered here covers a wide variety of stellar masses and 
luminosit ies, as in the observe d one. Nevertheless, as dis- 
cussed bv lGonzalez et alj (|2009T ), Shankar et al. (2009c) and 
further below, the model has a tendency to produce broader 
distributions at fixed size and/or stellar mass than those ac- 
tually observed. We believe that this effect is mainly due to 
somewhat inappropriate physical recipes more than inade- 
qu ate selection cuts in t he model (see the detailed analysis 
bv lGonzalez et al"1l2009l on this exact point). 

We identify the formation epoch of an early-type galaxy 
as the first time along the merger tree that the progenitor 
becomes an early-type. Note that this classification of early- 
types and their formation does not make any assumptions 
about the color or the star formation rate. In the following, 
we will present results based on the full sample of early-type 
galaxies, regardless of whether they end up being central 



or satellite galaxies at z = 0. However, because we have 
removed BCGs from the data, we have checked that our 
basic result, of a flat size-age relation at fixed stellar mass, 
is still conserved if we remove the objects which are centrals 
at z = in the model halos. 

3.1 Higher densities at early times? 

Before we present a more direct comparison of these models 
and the data shown in the previous figure, it is interesting to 
see if the models are consistent with the notion that objects 
at high redshift are denser, for the reasons discussed in the 
Introduction. Figure shows the sizes (top) and densities 
(bottom) as a function of lookback time, for objects of fixed 
stellar mass Mstar at each lookback time. The three panels 
show these predictions for three choices of Mstar. Notice 
that high-redshift objects are smaller (and hence denser) 
than their counterparts of the same Mstar at later times - 
in qualitative agreement with observations and expectations. 

However, notice that the models predict the typical size 
to decrease as Mstar increases. This is grossly discrepant 
with observations, and suggests that there is considerable 
room for improvement with rega rds to how the models a ssign 
sizes - this was also noted by IGonzalez et al.l (|2009h and 
Shankar et al. (2009c). We further discuss in Shankar et 
al. (2009c) the actual successes and failures of the present 
model by making use of a combined comparison with the size 
and mass distributions as observed in SDSS. In this paper, 
however, we are less interested in the form of the correlation 
between size and mass — we are more interested in checking 
if a "chaotic" formation model, such as a hierarchical one, in 
which early-type galaxies form and grow continuously, can 
reproduce the fact that the size-mass relation is independent 
of age - as observed for today's massive galaxies, irrespective 
of whether the final size is actually the one observed for that 
bin of stellar mass. 
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3.2 Mass dependent evolution of sizes and masses 

Figure [6] shows the model predictions in a format that is 
closer to that shown in the previous section. We here plot the 
median sizes, stellar masses and densities (along with their 
dispersions), averaged over 100 merger historieo of galaxies 

We have checked that the main results of this paper do not de- 
pend significantly on the specific number of random merger his- 
tories adopted. Increasing this number by a factor of a few yields 



with stellar mass at z = in the mass bin indicated at the 
top of each column. In each case the trees are followed back 
in time, choosing the most massive early-type progenitor, 
until this is no longer possible. The blue squares refer to the 
properties of galaxies (size, mass, or density) at z = 0, while 
the red circles refer to the properties the same set of galaxies 



essentially equal results both in the mean trends and broadness 
of distributions around the mean. 
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had at their formation epoch, if orm (defined as lookback time 
in Gyr). The blue squares in the top panels, which show 
the z — size-age relation for galaxies in narrow bins in 
Mstar, can be directly compared with the measurements 
shown in Figure [3] assuming the age of a galaxy is a good 
proxy for its formation epoch tf or m- Notice that, at fixed 
Mstar, the size-age relation is weak, with a slight tendency 
for old objects to be smaller in the smallest Mstar bin, but 
to be larger in the largest Mstar bin. Indeed, the models 
suggest that, at small A/star, older objects have slightly 
smaller sizes whereas the opposite is true at higher Mstar- 
The sense of these weak trends at fixed Mstar is also in 
qualitative agreement with our measurements in the SDSS, 
and may be considered a significant success of the models. 
However, we stress that the models predict smaller sizes at 
large Mstar, whereas the data shows the opposite trend. As 
a result, the overall R e — Mstar relation in these models is 
grossly discrepant with that in the SDSS. (This is essentially 
the same problem we found in Figure [5] the only difference 
is that there the galaxies in a given panel were selected to 
have the same Mstar at all lookback times, whereas here 
they have the same Mstar only at z = 0, and we then study 
their progenitors at earlier times.) 

If the problem with the models is simply the overall 
normalization of the sizes for a given Mstar, then it is in- 
teresting to study what aspect of the models produced the 
flat size-age relation. The first question which arises in this 
context is what these objects looked like in the past. The 
red circles in each panel show the sizes, masses and densi- 
ties (top to bottom) of these objects at the time they formed. 
The middle panels show that the mass change is larger for 
objects which formed at larger redshift, as one might expect, 
but that this change is most dramatic for the oldest most 
massive galaxies (which have increased their mass by about 
a factor of ten). Comparison with the top panels shows that 
the sizes increase when the masses do, and that, except in 
the low mass bin, the fractional increase in size since forma- 
tion is rarely larger than that in the mass. As a result, the 
densities of low mass galaxies today have decreased since 
they formed, whereas this trend is less clear for the higher 
mass galaxies. 

Figure [7] compares the models and the data in a format 
which is more like Figure [3] size versus formation time for a 
few bins in stellar mass (we again here make the assumption 
that the age of SDSS galaxies is a good proxy for their forma- 
tion time iform). To circumvent the problem that the scaling 
of SAM size with stellar mass is wrong, we show the sizes 
in a fixed mass bin normalized by the mean sizes of galaxies 
which formed most recently. For the SDSS data (open dia- 
monds), this means that we take the values shown by the 
diamonds in Figure[3]and divide by the value of the left-most 
diamond (for the galaxies with mass above 3 x 1O 1O M , we 
consider only galaxies with ages above ~ 3 — 4 Gyr, given 
that the bins below are not statistically significant); for the 
SAMs (filled circles), we take the blue circles shown in Fig- 
ure [6] divided by the leftmost blue circle (for consistency, for 
the more massive galaxies in the SAM we only select objects 
with a minimum age of ~ 3 — 4 Gyr) . Note how the SAMs 
are generally in good agreement with the data. They pro- 
duce rather flat size-age relations, and, at low masses, the 
older galaxies are about 2x smaller than younger ones. We 



show below that, in the models, this happens because lower 
mass galaxies undergo fewer mergers. 

3.3 Major versus minor mergers 

Figure [8] shows the predicted mean number Amergers of 
wet (dotted) and dry (solid), minor (left) or major (right) 
mergers that today's early-type galaxies have undergone 
since they formed. More specifically, we compute the mean 
number of mergers per Gyr a galaxy had since its formation 
epoch, averaged over all galaxies, as a function of lookback 
time t. (The numbers of mergers for each galaxy were ex- 
tracted from the same 100 merger trees used in the previous 
figures; wet mergers have a cold-total gas mass fraction in 
the progenitors that is greater than 0.15, else the merger 
is dry; minor mergers have mass ratios < 1 : 3, else the 
merger is major.) Notice that there are essentially no wet 
mergers; massive objects have had an order of magnitude 
more major mergers and at least a factor of two more minor 
mergers, than lower mass objects; and that minor mergers 
are typically a factor of 5 times more frequent than major 
mergers. 

Table 1 summarizes the actual increases in size and stel- 
lar mass experienced by the early-type galaxies in the model. 
It is apparent that while a substantial fraction of the stellar 
mass is added to the galaxy via major mergers, the sizes 
mostly increase via minor dry mergers since their forma- 
tion epoch. These galaxies remain gas-poor for most of their 
assembly history: only < 10% (decreasing to < 4% for the 
most massive systems) of the final sizes and masses grow due 
to gas-rich mergers. These objects are all bulge-dominated, 
given that < 10% of the size and stellar mass is increased 
via disk instability. 

The minor merger-dominated size evolution of massive 
early-type galaxies in Fig ure[8]may be related to that of their 
DM halos. lStewart et all (|2008l ) from high-resolution ACDM 
N-body simulations found that the mass assembly in "galac- 
tic" halos, those with mass in the range 10 n -10 13 /i _1 M , 
to be dominated by mergers that are ~ 10% of the fi- 
nal halo mass. Minor dry mergers would also more eas- 
ily preserve the proje ctions of the fundamental plane, such 
as th e the L-a (e.g., [Fabcr & J ackson! Il97fj; iDavies et all 
I l983h and Rr-L (e.g .. lKormendvlll977l ; IZiegler et all J 19991 ; 
Bernardi et alj |2003') relations (see discussion in lCiottill2009l ; 
lBernardill20091 and references therein) . Preliminary measure- 
ments also show that the central densities within the same 
physical scale for samples of low and high redshift galaxies of 
the s ame stellar mass, are consistent within a f actor of ~ 2 
Ce.g.. lCimatti et al.|[2008l ; IBezanson et alj|2009h . The latter 
findings could be consistent with an inside-out evolutionary 
scenario, where stellar matter is continuously added to the 
outskirts of the compact high-redshift galaxies as time goes 
on. However, larger samples of galaxies at different redshifts 
with well measured density and metallicity radial profiles 
are r e quired to set definite conclusions (e.g.. ICimatti et al.l 
l200St IBezanson etafl 120091; iHopkins et all l2009al ). In par- 
ticulai~ TBezanson et al.l (120091 . see also Shankar & Bernardi 
2009) discussed the results of several basic models for the 
size and mass evolution of spheroids. Overall, they conclude 
that galaxies with stellar mass ~ lO n M at z ~ 2, should 
undergo about 8 minor mergers to efficiently increase their 
sizes by a factor of ~ 5 and mass of ~ 2 to grow onto the local 
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Figure 8. Comparison between the mean number of wet and dry mergers per Gyr, averaged over 100 galaxies, extracted from the merger 
trees of the Bower et al. (2006; blue lines) catalogues. Each row shows the merger history, averaged over 100 realizations, of galaxies with 
stellar mass at 2 = in three different mass bins, as labeled. In the left column we plot the mean number of minor mergers, with mass 
ratio < 1 : 3, while the right column shows the mean number of major mergers with mass ratio > 1 : 3. The dotted and solid lines refer 
to the mean number of wet and dry mergers, defined to have a (cold) gas-to-total mass fraction in the progenitors higher and lower than 
0.15, respectively. 



Table 1. Fractional Increase In Stellar Mass and Radius. 





10 < logMstar < 10.5 


10.5 < logM star < 11 


11 < lc 


)gM st ar < H.5 


Type 


R MsTAR. 


R Mg TAR 


R 


MsTAR 



Minor Dry Mergers 57% 33% 73% 30% 90% 36% 

Major Dry Mergers 26% 49% 22.5% 59% 0% 61% 

Wet Growth 10% 8% 4% 2% 2.5% 1% 

Disk Instability 7% 10% 0.5% 9% 7.5% 2% 



Notes: Fractional increase in size and stellar mass for the galaxies of different mass at z = 0; although most of the stellar mass is added 
via major dry mergers, most of the growth in size is through minor dry mergers. 



size-mass relation. They also note that their central densities 
would then be consistent (within a factor of ~ 2) with the 
ones of SDSS galaxies with mass a factor of ~ 2 higher, fur- 
ther supporting the minor merger hypothesis. However, the 
cumulative number of minor mergers in the B06 model is sig- 
nifican tly lower than the one put forward bv lBezanson et al] 
{200$ ) (see also Shankar et al. (2009c)), and in fact the 
size increase in the mock massive galaxies presented here 
is hardly enough to bring them onto t he local size-mass re - 
lation, as discussed above. Moreover, iNipoti et al~l (|2009T l 
recently concluded through a series of N-body simulations, 



that reproducing the growth of a factor of about ~ 5 only 
through dry mergers is problematic. They find mergers not 
to efficiently grow galaxies in the required proportion, and 
to increase the scatter in the galactic scaling relations be- 
yond what allowed by observations. They therefore conclude 
that observational biases in the measurement of the com- 
pact high-z galaxy sizes, coupled to extreme fine-tuning in 
the merger processes, are required to accommodate a pure 
merger-driven scenario as main driver for the size evolution 
of ellipticals. 
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4 DISCUSSION 

While we were completing this wo rk, we became aware of th e 
study of the age-size relation by Ivan der Wei et all [|2009h . 
They too find that SDSS galaxies show no relation between 
size and age. Using simple prescriptions for the merger his- 
tories of galaxies between their formation redshift and the 
present, they also conclude that models in which galaxies 
grow through dry mergers are consistent with the observed 
evolution since z ~ 2 in the mean size and in the co- 
m oving mass density. H owever, as also recently addressed 
bv lBernardi et all ([2009), there is one important respect in 
which our results differ from theirs. 

Although we have studied the size as a function of age, 
their Figure 1 shows the age as a function of size. To en- 
able comparison with their work, our Figure [5] shows this 
relation at fixed Md yn and a in our dataset. While we agree 
with them that, at fixed a, this relation is weak, we come 
to a somewhat different conclusion about this relation at 
fixed Mdyn- Whereas they find that, at fixed Md yn , smaller 
galaxies are older, we find no correlation between age and 
size at fixed Md yn . This is consistent with our other results 
above, if one allows for the fact that Mstar and Md yn arc 
closely related, so it is a reasonable appr o ximat ion to sub- 
stitute one for the other. iBernardi et al.l (|2009l ) show that 
this discrepancy is almost certainly due to the fact that our 
sample is less contaminated by objects with disks: whereas 
ellipticals have a flat age-size relation at fixed Mdyn, age 
and size are anti-correlated for SOs and Sas. This suggests 
that the two early-type galaxy populations have different 
formation histories. 
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5 CONCLUSIONS 

We studied the size-mass-age relations as derived from a 
sample of about 45,700 early-type galaxies selected from 
SDSS, and compared our results with Bower et al.'s (2006) 
model of hierarchical galaxy formation. Our results can be 
summarized as follows. 

• At stellar masses below 3 x 1O 1O M0, the effective ra- 
dius i? e can be up to a factor of ~ 2 lower for older galaxies. 
However, more massive galaxies all share similar size distri- 
butions, irrespective of their ages (Figure[3)|. At these higher 
masses, the scatter in sizes at fixed stellar mass is about a 
factor of 4 (~ 0.6 dex), decreasing with increasing stellar 
mass (Figure [1}. These findings are at variance with a pure 
passive evolution, which would predict older galaxies to be 
much more compact at fixed stellar mass. 

• Hierarchical galaxy formation models predict that 
galaxies of the same stellar mass formed at high redshifts 
are much more compact and dense (Figure [5]), in agreement 
with observations at higher redshift (see §[T]). This is both 
because the high redshift Universe itself is denser, and be- 
cause dissipation effects are more effective at early times. 

• SAMs based on a hierarchical growth of galaxies driven 
by a first major, wet merger and a sequence of late, mi- 
nor, dry mergers, predict that these extremely small, high- 
redshift galaxies can grow, on average, onto the same lo- 
cal size-age relation (Figures [5] and [7]) . Note, however, that 
(dry) mergers are not the unique way to inc rease early-type 
galaxy sizes. For example. iFan et al.l (|2008t ) . put forward a 



Figure 9. Age-size relations at fixed dynamical mass Md yn (up- 
per panels), and at fixed velocity dispersion a (lower panels). 



model that postulates a strong galaxy expansion caused by 
the mass loss due to quasar feedback and stellar winds. This 
model predicts a local size-age relation that is consistent 
with that one observ ed in SDSS, at least at lowe r masses (see 
further discussion in IShankar fc Bernardjl2009h - something 
which the SAMs are unable to accomplish. 

• In the SAMs, galaxies which form at higher redshifts 
experience more mergers than galaxies which formed more 
recently, increasing their original sizes by a greater factor 
than galaxies which formed later. This process almost com- 
pletely wipes out the monolithic effect, growing all galaxies 
towards the same size-mass relation today. In particular, 
although most of the stellar mass is added via major dry 
mergers, most of the growth in size is through minor dry 
mergers. Minor mergers (mass ratios < I : 3) outnumber 
major mergers by about a factor of 5 at Mstar> I0 10 ' 5 Mq, 
and by a factor of 10 at smaller masses (Figure [5] and Ta- 
ble [TJ. 

• However, these SAMs provide a poor match to the 
local size-mass relation, and much more work has to be 
done to understand the origin of thes e discrepancies (e.g. 
iTruiillo et ail 120091 . iTavlor et all 120091 . and Shankar et al. 
2009c). 
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